A Topology— Changing First Order Phase Transition 



and the 



Dynamics of Flavour 



Tameem Albash, Veselin Filev, Clifford V. Johnson^? Arnab Kundu 



Department of Physics and Astronomy 
University of Southern California 
Los Angeles, CA 90089-0484, U.S.A. 



albash, filev, johnsonl, akundu Ousc.edu 



Abstract 

In studying the dynamics of large N c , SU(N C ) gauge theory at finite temperature 
with fundamental quark flavours in the quenched approximation, we observe a first order 
phase transition. A quark condensate forms at finite quark mass, and the value of the 
condensate varies smoothly with the quark mass for generic regions in parameter space. 
At a particular value of the quark mass, there is a finite discontinuity in the condensate's 
vacuum expectation value, corresponding to a first order phase transition. We study the 
gauge theory via its string dual formulation using the AdS/CFT conjecture, the string 
dual being the near-horizon geometry of N c D3-branes at finite temperature, AdSs- 
SchwarzschildxS 15 , probed by a D7-brane. The D7-brane has topology K 4 x S 3 x S 1 and 
allowed solutions correspond to either the S 3 or the S 1 shrinking away in the interior 
of the geometry. The phase transition represents a jump between branches of solutions 
having these two distinct D-brane topologies. The transition also appears in the meson 
spectrum. 
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1 Introduction 



String theory is a powerful tool for probing the strongly coupled dynamics of gauge theory, 
physics which is of vital importance especially in the context of the strong nuclear interactions. 
Gauge/string correspondences have enlarged and refined the toolbox available for such studies, 
and there is a large literature on the subject, with several powerful examples such as the 
AdS/CFT correspondence [U [21 E] and deformations thereof [4|. 

However, we are still some way from describing the "realistic" dynamics of QCD using this 
approach. The main challenges that remain include getting access to low N c , fully including 
dynamical quarks in the fundamental representation of SU(N C ), and getting reliable control of 
the non-supersymmetric regime. 

It is to be hoped (if not expected) that even if we cannot obtain a controllable string dual of 
QCD, there may be considerable progress to be made in capturing physical phenomena that 
are in the same universality class as those of QCD. This is the motivation of the present work, 
which is part of a series of studies upon which we hope to report new and interesting results. 

We study the geometry of AdSs-Schwarzschild xS 5 , which is the decoupled/near-horizon ge- 
ometry of N c D3-branes, where N c is large and set by the (small, for reliability) curvature of the 
geometry. The physics of closed type IIB string theory in this background is dual to the physics 
of M = 4 supersymmetric SU(N C ) gauge theory in four dimensions, with the supersymmetry 
broken by being at finite temperature [5]. The temperature is set by the horizon radius of the 
Schwarzschild black hole, as we will recall below. 

We introduce a D7-brane probe into the background. Four of the brane's eight world-volume 
directions are parallel with those of the D3-branes, and three of them wrap an S 3 C S 5 . The 
remaining direction lies in the radial direction of the asymptotically AdSs geometry. 

Such a D3-D7 configuration controls the physics of the SU(N C ) gauge theory with a dynamical 
quark in the fundamental representation [Bj. The configuration (at zero temperature) preserves 
M = 2 supersymmetry in D = 4, and the quark is part of a hypermultiplet. Generically, we 
will be studying the physics at finite temperature, so supersymmetry will play no explicit role 
here. 

We are studying the D7-brane as a probe only, corresponding to taking the N c ^> Nf limit, 
and therefore there is no backreaction on the background geometry. This is roughly analogous 
to the quenched approximation in lattice QCD. The quark mass and other flavour physics-such 
as the vacuum expectation value (vev) of a condensate and the spectrum of mesons that can 
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be constructed from the quarks-are all physics which are therefore invisible in the background 
geometry. We will learn nothing new from the background; our study is of the response of the 
probe D7-branes to the background, and this is where the new physics emerges from. 

We carefully study the physics of the probe itself as it moves in the background geometry. The 
coordinates of the probe in the background are fields in an effective D7-brane world-volume 
theory, and the geometry of the background enters as couplings controlling the dynamics of 
those fields. One such coupling in the effective model represents the local separation, L(u), of 
the D7-brane probe from the D3-branes, where u is the radial AdSs-Schwarzschild coordinate. 

In fact, the asymptotic value of the separation between the D3-branes and D7-brane for large 
u yields the bare quark mass m q and the condensate vacuum expectation value (vev) (ipip) as 
follows [3 [S]: 

lim Liu) = m H — - + ... (1) 

where m = 2iTa'm q and — c = (ipilj)/(87r 3 a'NfTf), where in this paper Nf = 1. (The fun- 
damental string tension is defined as T = l/(27ra') here, and the D7-brane tension is r 7 = 
(27r)~ 7 (o/) -4 .) The zero temperature behaviour of the D7-branes in the geometry is simple. 
The D7-brane worldvolume actually vanishes at finite u, corresponding to the part of the brane 
wrapped on the S 3 shrinking to zero size. The location in u where this vanishing happens en- 
codes the mass of the quark, or equivalently, the separation of the probe from the D3-branes. 
In addition, in the zero temperature background, the only value of c allowed is zero, meaning 
no condensate is allowed to form, as is expected from supersymmetry. 

The finite temperature physics introduces an important new feature. As is standard[9J, finite 
temperature is studied by Euclideanizing the geometry and identifying the temperature with 
the period of the time coordinate. The horizon of the background geometry is the place where 
that S 1 shrinks to zero size. The D7-brane is also wrapped on this S l , so it can vanish at the 
horizon, if it has not vanished due to the shrinking of the S 3 . For large quark mass compared 
to the temperature (horizon size), the S 3 shrinking will occur at some finite u > Uh, and the 
physics will be similar to the zero temperature situation. However, for small quark mass, the 
world-volume will vanish due to the shrinking of the S l corresponding to the D7-branes going 
into the horizon. This is new physics of the flavour sector. 

The authors of ref.[10j explored some of the physics of this situation (the dependence of the 
condensate and of the meson mass on the bare quark mass), and predicted that a phase transi- 
tion should occur when the topology of the probe D7-brane changes. However, they were not 
able to explicitly see this transition because of poor data resolution in the transition region, 
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coming from using UV boundary conditions on the scalar fields on the D7-brane world-volume. 
The origin of this phase transition, as we shall see, is as follows: The generic behaviour of an 
allowed solution for L(u) as in equation ([I]), is not enough to determine whether the behaviour 
corresponds to an S" 3 -vanishing D7-brane or an S^-vanishing D7-brane. The choices of branch 
of solutions have different values of c, generically. In other words, for a given value of m there 
can be more than one value of c. There are therefore two or more candidate solutions po- 
tentially controlling the physics. The actual physical solution is the one which has the lowest 
value for the D7-brane's free energy. The key point is that, at a certain value of the mass, the 
lowest energy solution may suddenly come from a different branch, and, as the corresponding 
value of the condensate changes discontinuously in moving between branches, we find that the 
system therefore undergoes a first order phase transition. On the gauge theory side, we can 
imagine a similar situation occurring; two different branches of solution are competing, and the 
lowest energy branch is always picked. We are able to uncover this physics by doing a careful 
numerical analysis of the equations of motion for the probe dynamics on the gravity side of 
the AdS/CFT correspondence, which goes well beyond that carried out in previous papers. 
By using IR boundary conditions instead of UV boundary conditions, our analysis allows us 



to obtain vastly more data points, allowing us to complete the picture described above 3- In 
addition, we are able to study the meson spectrum in some detail, and we find that the phase 
transition also manifests itself as a first order discontinuity in that physics. 

2 The Probe Computation 

We begin by reviewing the physics of the D7-brane probe in the AdSs-Schwarzschild back- 
ground solution[T0]. The metric is given by: 

d s 2 = -^-dt 2 + -^--du 2 + ^-dx- dx + R 2 d6 2 + R 2 cos 2 9dQ 2 3 + R 2 sin 2 0# 2 , (2) 
R l j (u) R 2 



where x is a three vector, 



and the quantity R is given by: 



b 4 

/(«) = u 2 -- 2 , 
u 



R 2 = ^47rg s N c a' 



^dn the process of writing this paper, similar work was published by Mateos et al.[ll]. In addition, after this 
paper appeared, I. Kirsch pointed out to us some related work in his thesis [12]. 
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where g s is the string coupling (which, with the inverse string tension a' sets for example, 
Newton's constant). The quantity b is related to the mass of the black-hole, b 2 = 8G5mb.h./(37r). 
The temperature of the black hole can be extracted using the standard Euclidean continuation 
and requiring regularity at the horizon. Doing this in the metric given by equation (J2J), we find 
that = b/%R 2 . Therefore, by picking the value of b, we are choosing at which temperature 
we are holding the theory. We choose to embed the D7-brane probe transverse to 9 and <fi. 
In order to study the embeddings with the lowest value of the on-shell action (and hence the 
lowest free energy), we choose an ansatz of the form = and 9 = 9(u). The asymptotic 
separation of the D3 and D7-branes is given by L(u) = usm9. Given this particular choice of 
embedding, the worldvolume of the D7 brane is given by: 

^T g = u 2 cos 3 9{u) Vdet S V« 2 + (u 4 - b A ) 9'{u) 2 , (3) 

where g is the determinant of the induced metric on the D7-brane given by the pull-back of 
the space-time metric G^. We are interested in two particular cases. First, there is the case 
where u goes to b, which corresponds to the the D7-brane probe falling into the event horizon. 
In the Euclidean section, this case corresponds to the shrinking of the S l of periodic time. 
This corresponds to what in ref.pl)] were called condensate solutions. Second, there is the 
case where 9 goes to tt/2, which corresponds to the shrinking of the S 3 . This corresponds to 
what were called Karch-Katz-like solutions in ref. [10]. It is this change in topology (S 1 versus 
S 3 shrinking) between the different solutions that will correspond to a phase transition. The 
classical equation of motion for 9{u) is: 

d f u 2 (w 4 — & 4 ) 9'{u) cos 3 (9{u)) \ 9 „„. . „. , /-^ — , AS „ , Nn , \ 

Tu[ ^ + y _ W) m J + 3 co (M) 3in { W + ( ' m = °' (4) 

When u goes to infinity and the background metric becomes asymptotically AdS 5 x S 5 , the 
equation of motion reduces to: 

4- (u 5 9'(u)) + 3u 2 9(u) = 0, (5) 
du 

which has solution: 

9(u) = -fm+4) ■ (6) 

These two terms are exactly the non-normalizable and normalizable terms corresponding to 
a dimension 3 operator (ipip) in the dual field theory with source m and vacuum expectation 
value (vev) c. 
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3 Embedding Solutions 



We solve equation (j3J) numerically using a shooting technique. Shooting from infinity towards 
the horizon, physical solutions are those that have a finite value at the horizon. This will 
only be accomplished for a particular m and c value from equation (JTJ) , which would have to be 
delicately chosen by hand. Therefore, if we instead start from the horizon with a finite solution, 
it will shoot towards the physical asymptotic solutions we desire. In order to be able to analyze 
the phase transition between the condensate and Karch-Katz-like solutions, we shoot from 
the horizon for the condensate solutions and from 9 = 71/ 2 for the Karch-Katz-like solutions. 
This technique avoids having to correctly choose the boundary conditions at infinity, which is a 
sensitive procedure, allowing us to have many more data points to analyze the phase transition. 
We impose the boundary condition that, at our starting point-the horizon, we have: 

Wls^o = -g-taa0(6) 

e'(u)\ s ^ = 00 (7) 

We argue that this is the physical boundary condition to take; the first is simply a result of 
taking the limit of u — > b in the equation of motion in equation (HI) , whereas the second is a 
result of requiring no conifold singularity as the S 3 shrinks to zero size 



We solve the equation of motion, equation (jlj), numerically using 6 = 1 and R — 1. These 
numerical choices correspond to fixing the temperature and to measuring lengths in units of 
the radius of the AdS space; the latter condition also means we are choosing a particular 
relationship between the t'Hooft coupling and the dual quantities: 

^ = 9 2 ymN c = ^- 2 , /r 1 ^. (8) 

Several D7-brane embedding solutions are shown in figure [U The red (solid) lines correspond 
to Karch-Katz-like solutions, and the blue (dashed) lines correspond to condensate solutions. 
From each of these solutions, we can extrapolate the bare quark mass and quark condensate 
vev. 



4 The Phase Transition 



We plot the c values as a function of m in figure 2(a) When enlarged, as shown in figure 2(b) 



we find, as anticipated in the introductory remarks, the multi-valuedness in c for a given m. 
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Figure 1: Solutions for the D7-brane probe in the black hole background. 



Physics will choose just one answer for c. There is therefore the possibility of a transition from 
one branch to another as one changes m. 

In order to determine exactly where the transition takes place, we have to calculate the free 
energy of the D7-brane. In the semi-classical limit that we are considering, the free energy is 
given by the the on-shell action times f3~ l . For our case, this is simply given by [8]: 



P 



^r 7 Nf J d A x dil 3 du a/— det g , 



(9) 



where here Nf = 1. We calculate this integral numerically using our solutions for 6(u), and 



we plot the results for the energy in figure 2(c) after we regulate the result by subtracting off 
the energy from the 8 = solution. There is again multivaluedness at the same value of m as 



before, and we zoom in on this neighbourhood in figure 2(d) If one follows the solutions of 
lowest energy for a given m, one can clearly see that there is a crossover from one branch to 
another as m changes. This was also observed in ref. [14J. Therefore, we find a first order phase 
transition — at m « 0.92345 — where the condensate's vev jumps discontinuously. This had 
been deduced on other grounds in ref. [15]. We show where the jump between curves occurs 



with the dashed green line in figures 2(b) and 2(d) . We note that the shaded areas in figure 2(b) 
are equal; this is to be expected since c and m are thermodynamically conjugate to each other, 
and the area of the graph has the interpretation as a free energy difference dS ~ —cdm. 
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5 Meson Spectrum 

Next we study the meson spectrum. This can be read off from the physics of the fluctuations in 6 
and about our classical solutions #(o) = 9{u) and 4>(o) = described above [16j, corresponding 
to scalar and pseudo scalar fields respectively in the gauge theory [TT]. As a reminder, in ref. [TT] . 
the exact meson spectrum for the AdSs x S* 5 background was found to be given by: 

M(n,£) = ^^( n + £+l)(n + £ + 2) , (10) 

where £ labels the order of the spherical harmonic expansion, and n is a positive integer that 
represents the order of the mode. In our meson spectrum calculations, we have only considered 
the n = £ = solutions for the fluctuations in 9 and </>. We begin by considering the four- 
dimensional mass of the meson corresponding to fluctuations in 0. By considering an ansatz of 
the form: 

4>(z,t) = + S4>(z,t) = f(z)e-^, (11) 
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where z = u~ 2 in terms of our previous coordinates, and by expanding to second order in 
5<f)(u,t), we can calculate an equation of motion for f(z). In particular, we find that near the 
event horizon, the equation of motion becomes: 



f(z)+ j^L + ^i m = o (12) 



This equation has solutions given by in-falling and out-going waves: 

f(z) = a(l-zb 2 ) iS *r +b(l-zb 2 y iJ %r . (13) 

This analysis of the quasinormal mode behavior has been shown in refs. |181 [T9l 120] . and the 
(correct) boundary condition of in-falling waves corresponds to the fact that fluctuations about 
condensate solutions do not allow for mesons with a discrete mass spectrum^. The mesons 
"melt" away into the plasma, from the gauge theory perspective. However, as was shown 
in ref.[2n], we can still find the mass of the meson (which we shall denote as M) for these 
fluctuations before they melt, in addition to the meson's lifetime (r) with the identification: 

Re[uj] = M , Jm[u] = (2r) _1 . (14) 

The fluctuations about the Karch-Katz-like solutions have a discrete mass spectrum; therefore 
the mesons are stable, and u will be purely real. To search for the condensate fluctuations, we 
do a convenient field redefinition: 

f(z) = y{z)(l-zb 2 y iB ^ , (15) 

such that our boundary condition is simply: 

y{b- 2 ) = 1 

y'ib- 2 ) ~ (l_^)-i| rf2 ^oc. (16) 
For the fluctuations about Karch-Katz-like solutions, we simply use the boundary conditions: 



(*)=t/2 ' 



/'WLw^/a = °> ( 17 ) 

where e is sufficiently small. (We used e ~ 10~ 2 here.) In order to find the correct value 
of io for both sets of fluctuations, we note that in the limit that z approaches zero, where the 



2 Refs.[l9| [20] appeared after the original version of this manuscript appeared, and were useful in refining 
this section's discussion. 



9 



background becomes asymptotically AdS 5 x S 5 , the field f(z) must only be comprised of a 
normalizable term, i.e., the coefficient in front of the non-normalizable term must be zero. We 
show the first level of the discrete four-dimensional meson mass dependence on the bare quark 



mass in figures 3(a) and 3(b) (zoomed). In figure 3(c) we show the behavior of the real and 
imaginary part of uj. 

For large bare quark mass, our solutions match those of the AdS 5 x S 5 case discussed in ref. [TT] 
(shown as the black lines with long dashes). We note that, for large bare quark mass, the meson 
mass is bounded from above by the AdSs x S 5 mass. This is to be expected since the coupling 
in the finite temperature gauge theory goes as: 



1 



9ym 



1 



-1/2 



5ym 



/3=oo 



In addition, we see that at zero bare quark mass, the meson has a finite mass. This is simply the 
mass gap of the three-dimensional gauge theory: In the Euclidean section, the four-dimensional 
gauge theory is compactified on an S l , whose radius is set by the energy scale of the theory. In 
this case, the energy scale is set by the bare quark mass, so at zero bare quark mass, the circle 
shrinks to zero size, making the gauge theory effectively three-dimensional. 
For fluctuations in the 9 direction, we can follow a similar analysis as for the fluctuations in 
the (p direction. We begin by searching for the 4-dimensional mass by taking an ansatz = 
and 9 = 9/q) + 59(z,t), where S9(z,t) has the form: 



S9(z,t) =f(z)e 



-iujt 



(19) 



The equation of motion shows the same behavior as the 0-meson, so we choose the same 
boundary condition as in equation ( fT6l) for the condensate solutions. For the Karch-Katz-like 
solutions, we use the boundary condition: 



/(*)| 
f{z)\ 



W=*/2 
W=w/2 



OO 



(20) 



where e is suitably small. (We chose e ~ 10 .) The spectrum is shown in figures 4(a) and 



4(b) In figure 4(b), the disconnect between the condensate and Karch-Katz-like solutions is 
due to there being tachyonic modes that extend beyond the zero meson mass [TH] , but we have 
not shown them in our plots. They occur well away from the physical regime. 
We may also consider fluctuations in the gauge field living on the D7-brane; such fluctuations 
would correspond to vector mesons in the gauge theory. In order to study their behavior, we 
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(a) (b) 



Re [u] 




(c) 

Figure 3: Four-dimensional meson mass M, as a function of fundamental quark mass m, from 
fluctuations in 0. Only the first level of the discrete spectrum is shown. The blue (dashed) 
curves in (a) and (b) correspond to the meson mass before the meson melts. 



must begin by including the field strength term (from the world-volume gauge field) to the 
DBI action. The action is given by: 

S = -T 7 J d 8 £ (- det (g ab + 2ira'F ab )f /2 + ^p-fr J P [C (4) ] A F A F , (21) 

where P [C(4)l is the pull-back of the 4-form potential sourced by the N c D3-branes. Our 
previous solutions assumed a solution A b = , so we will consider quadratic fluctuations about 
this solution. Expanding this action to second order in A&, we have: 



S 



T7 / 



det (gat,)) 



-1/2 



1 - -(27ra') 2 F ab F ba 



8 R A 



ij-Tklt 



(22) 



where the subscripts a, b run over the world-volume coordinates and i, j, k, I run over the world- 
volume coordinates transverse to the D3-brane. We have dropped the first order term since 
it will just give the classical equation of motion for A b . Since we are only interested in the 
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Figure 4: Four-dimensional meson mass M, as a function of fundamental quark mass m, from 
fluctuations in 9. Only the first level of the discrete spectrum is shown. The blue (dashed) 
curves in (a) and (b) correspond to the meson mass before the meson melts. 



components of the gauge field along the D3-brane worldvolume, the Wess-Zemino term will 
not contribute to the equations of motion that we are interested in. Therefore, the equation of 
motion for the quadratic fluctuations is given by: 

d a (V^F a ») = . (23) 

To compute the 4-dimensional mass of the vector meson, we consider an ansatz of the form: 

A, = V,(z)e-^ . (24) 

We want to consider the case where the condition k^N 1 = still holds, which requires us to set 
At = since we are working in the rest frame of the vector meson. The equations of motion 
for the three V^i's are identical, so we only need to solve for one of them. The equations of 
motion show the same behavior as the 0-mesons, so we follow the same procedure as before 
with the same boundary conditions as in equations (fT6l) and (TTT|) . The spectrum is shown in 
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figures 5(a) and 5(b), We see again that the mass is bounded from above by the AdSs x S 5 
mass (shown as the black lines with long dashes). 




(c) 

Figure 5: 4-dimensional meson mass M, as a function of fundamental quark mass m, from 
fluctuations in A&. Only the first level of the discrete spectrum is shown. The blue (dashed) 
curves in (a) and (b) correspond to the meson mass before the meson melts. 



6 Conclusion 

We expect that, while the details of this construction will not persist in a "realistic" QCD string 
dual, the phase transition itself represents strongly coupled dynamics that may well persist as 
part of the full story of the QCD phase diagram. The transition is exciting in itself, of course 
(particularly since its dual involves a change of topology in the D7-brane world-volume), but 
much further work is needed on several questions. For example, the robustness of the phase 
transition against 1/N C corrections would be interesting to study. It would also be interesting 
to see if the back-reaction can be included, allowing us to study Nf ~ iV c and follow the 
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phase structure to this regime. Although there are complications involving conical deficits in 
otherwise flat directions when dealing with the back-reacted geometry in the case of a D3-D7 
system, some progress has been made in studying the back-reacted geometry of other more 
stable Dp-D(p+4) systems [2H 1221 1231 121 1251 123 123 • However , man Y of these calculations 
have not considered black hole or other finite temperature configurations. These, and several 
other questions, are exciting matters for further study. 
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